Short Range Smectic and Long Range Nematic Order in the Pseudogap Phase of 

Cuprates 



R.S. Markiewicz/ J. Lorenzana,^ G. Seibold,^ and A. Bansil^ 

^ Physics Department, Northeastern University, Boston MA 02115, USA 
^ISC-CNR, Dipartimento di Fisica, Universita "La Sapienza", Piazzale Aldo Moro 5, Roma, Italy 
^Institut Fiir Physik, BTU Cottbus, PBox 101344, 03013 Cotthus, Germany 

Wc present a model for the combined nematic and 'smectic' or stripe-like orders seen in recent 
scanning tunneling microscopy (STM) experiments in cuprates. We model the stripe order as an 
electronic charge density wave with associated Peierls distortion - a 'Pomeranchuk wave'. Disorder 
restricts this primary order to nanoscale domains, while secondary coupling to strain generates 
nematic order with considerably longer range. 

PACS numbers: 



Two of the most intriguing phases in strongly corre- 
lated matter are the Mott and high- Tc superconducting 
phases found in close proximity in cuprates. Any di- 
rect connection between these phases is blocked by the 
universal intervention of a pseudogap phase in the un- 
derdoped regime. Controversy swirls about this phase 
- is it associated with preformed pairs or a competing 
order, and if the latter, what is the nature of this 'hid- 
den order'? While there is considerable evidence for 
'smectic' or stripe-like order in underdoped cuprates [1, 2], 
such stripes should be sensitive to the quenched disor- 
der that is generally present in cuprates. Hence, long 
range stripe order is less likely than a frozen, glassy 
phase with only short-range orientational order. Alter- 
natively, there could be a fluctuating stripe phase con- 
sistent with a nematic order, which retains orientational 
symmetry-breaking while translational symmetry is re- 
stored. Totally unexpected therefore was the recent dis- 
covery in Bi2Sr2CaCu208 (Bi2212) of a coexisting short 
range (glassy) smectic order and a long range nematic 
order [3]. In Ref. 4, it was shown that fluctuations of 
the nematic phase are strongly correlated with the smec- 
tic phase, and this coupling was described in terms of a 
Ginzburg-Landau theory originally developed for liquid 
crystal phases. 

However, the origin of the nematic phase itself was not 
explained. In this Letter we provide an explanation for 
this smcctic-nematic phase in terms of a charge-density 
wave (CDW) coupled to strain. We assume that the pseu- 
dogap is associated with a competing density wave (DW) 
order, and that this order has been observed by scanning 
tunneling microscopy (STM) studies[5-12] of cuprates. 
Then the problem is tackled in two steps. First, we ana- 
lyze the occurence of charge order (CO) from the inter- 
play between electrons and lattice deformations within a 
microscopic model. In a second, more phenomenological 
approach, we incoporate the effect of impurities which 
break up the regular CO into small domains but still al- 
low for large scale nematic distortions in agreement with 
Ref. [3] (cf. Fig. 2). 

In mean-field type approximations of Hubbard mod- 



els the dominant stripe order appears to be predomi- 
nantly magnetic[13]. On the other hand exact diago- 
nalization of small clusters within the tJ-model supple- 
mented with electron-phonon coupling supports the view 
of stripes which are stabilized by lattice modes [14]. In 
fact, experimentally DW order is found to be accompa- 
nied by lattice distortions, at least in La2_a;Srj,Cu04_|_5 
(LSCO). Hence we here explore a model with a primary 
charge order, stabilized by a lattice distortion. Our con- 
siderations are based on the time-dependent Gutzwiller 
approximation applied to the one-band Hubbard model 
where the electron-lattice coupling is incorporated via a 
modulation of the hopping parameters. The calculations 
are summarized in the supplementary materials. Our re- 
sults are generic for any model of charge-ordered stripes 
coupled to a lattice distortion in particular, the main 
conclusions should also be valid for magnetically origi- 
nated stripes with induced charge order or any model 
of charge-ordered stripes coupled to strain. The resiilt- 
ing stripes provide a good description of many proper- 
ties of the modulated phase [5 11] found in STM stud- 
ies of Ca2Cu02Cl2 (CCOC), Bi2212, and BiaSraCuOe 
(Bi2201). 

We find that the stripes are stabilized by nesting 

of the flat, nearly parallel antinodal (AN) FS sections 
near (7r,0) in the Brillouin zone, [9-11, 15] and hence 
we refer to them as AN nesting (ANN) stripes. This 
ANN phase breaks C4 symmetry and hence also repre- 
sents a quasi-one-dimensional charge order (CO) state. [8] 
Figure 1 shows that the model correctly describes the 
mcasured[9, 10] doping dependence of the ANN peri- 
odicity. Figure 1 plots the calculated doping depen- 
dence of the ANN charge nesting vectors for the Bi 
cuprates, [16, 17] displaying a strong material depen- 
dence. Shown also are the experimental superlattice pe- 
riodicities for charge order in Bi2201[10] (green circles) 
and-2212[9] (violet circles). Clearly the experimental su- 
perlattices in the Bi-cuprates are close to the predicted 
ANN periodicities. However, it should be noted that ex- 
perimental data arc extracted at the energy scale of the 
pseudogap whereas the nesting vectors signal the transi- 
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FIG. 1: (Color online.) Nesting vectors for ANN (q, q) charge 
order, for Bi2201 (green solid line), Bi2212 with (violet dot- 
dashed) or without (violet dotted) bilayer splitting, compared 
to experimental (0, q) charge stripe vectors for Bi2201 (dark 
green circles) [10] and Bi2212 (violet circles) [9]. Band param- 
eters for the theoretical calculations are from Refs. 16 and 17 
respectively. 



tion towards CO at zero frequency. This might be the 
reason why the calculated ANN occurs along the diag- 
onals in contrast to the electronic order along the Cu- 
O bond direction oberved in STM. In Bi2212 there are 
two nesting vectors, associated with bonding and anti- 
bonding combinations of the bilayer-split bands, and the 
experimental data fall close to the bonding band nest- 
ing vector. In decply-underdoped Bi2212, this charge 
density wave (CDW) may be unstable against nanoscale 
phase separation. [18] Hints for this latter effect may have 
been observed in recent STM studies: Ref. 11 found that 
the phase we here identify as ANN seems to weaken be- 
low l/8th doping, while in a similar doping range in 
CCOC Ref. 12 found that islands of a phase without C4 
symmetry breaking became more prevalent with reduced 
doping. [19, 20] 

Allowing for the symmetry-broken CO phase it is found 
that the model can naturally explain a number of fea- 
tures of the pseudogap regime seen in STM and ARPES 
experiments. This is exemplified for an ordered stripe 
array (periodicity 6 lattice constants) where we refer to 
the Supplementary material for details of the correspond- 
ing analysis. In particular, we find that (1) there is 
a contrast reversal of the stripes in STM images, be- 
tween positive and negative energies, consistent with 
experiment[3, 8, 9, 21], Figs. S3(a,b); (2) maps of the 
spectral intensity at the Fermi energy typically show 
arc-like Fermi surfaces (FSs), consistent with ARPES 
results, which lose intensity close to the conventional 
(tt, 7r)-antiferromagnetic (AF) zone boundary as seen in 
STM[9], Figs. S3(f,g); (3) by including a secondary mag- 
netic order parameter, it is possible to get strongly angu- 
lar FS sections, which had been seen in ARPES and taken 
as evidence for underlying stripe phases[22]. Fig. S3(g), 



or even half-pockets^ which close on the same side of 
the AF zone boundary. Fig. S3(e), as seen in ARPES 
studies[23-26], although (7r/2, 7r/2) centered pockets may 
have a structural origin associated with low-temperature 
tetragonal order[27], unrelated to stripe physics; (4) 
ARPES finds that the pseudogap in the antinodal nesting 
regime is centered at the Fermi level [28], and this arises 
naturally in the ANN model, since it is precisely the AN 
FS sections that control the nesting. Fig. S4. 

However, the most puzzling feature of the pseudogap 
phase is the coexistence of a stripe- like order with a g = 
order, which breaks the symmetry of a unit cell - most 
likely by making the oxygens inequivalent.[3, 29] Such an 
extra symmetry breaking arises naturally in any verti- 
cal stripes, which have no bond length modulation along 
y (the direction along the stripe). Since the hopping 
parameter t scales as a"^, where a is the lattice con- 
stant and p ~ 3.5 for nearest neighbors, then when a 
is modulated by 5a the linear corrections average out 
over the stripes, but there is a finite quadratic correc- 
tion, 5tx ^ {5aY . Hence the Pomeranchuk wave couples 
quadratically to a 5 = shear mode, t^ > ty.[30] This in 
turn leads to a nematic electronic order: the hole doping 
is inequivalent on the in-plane x and y oxygens. [31, 32] 
This doping asymetry was recently found [3] in STM stud- 
ies of Bi2212 and is also consistent with the g = mag- 
netic order found in the pseudogap state of cuprates,[29] 
as long as the oxygen holes have a magnetic moment. 
Furthermore, if the stripes fluctuate in time or space (due 
to impurity pinning), then 5a = 5aocos{qxx) will average 
to zero, but since 5a^ = 5aQ[cos{2qxx) + l]/2 the shear 
contribution at q — will not. 

Within our model, nematic order is a true emergent 
phenomenon. The primary instability is to a CDW or- 
der, but in two dimensions this is unstable in the pres- 
ence of charged impurities [33]. On the other hand, there 
is a secondary coupling to a strain distortion, and this 
being long-range is much less sensitive to impurities. [34] 
Thus stripe domains will be disordered with x and y do- 
mains, while shear domains will be more robust, having 
long-range interactions, so the nematic domains will have 
larger correlation lengths, as seen in experiment. Here we 
provide a 'proof of principle' that the secondary order pa- 
rameter can have a larger correlation length. The most 
general Landau-Ginzburg effective Hamiltonian to model 
the strain-density wave interaction is: 

^eff ^ ^<P + + Hri-tj,, (1) 

where the DW order is characterized by a charge density 
modulation as [35] 

p(r) - p + [<^,(r)e*«==- + 0^(r)e'««^ + c.c], (2) 

and (px, 4>y are the two competing density waves in or- 
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thogonal directions. Then 



\d.^y\' + \dyq^.ni3) 



The strain Hamiltonian is[36, 37] 
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where the strain components are rj = (e^i 



-yy 



(5) 

ei = {exx + Eyy)/V^i 62 = Cxy [Since this is a one-band 
model, the effect of the strain on the oxygens, leading to 
nematic order, is implicit.] 

To compare with STM results, we can simplify the 
above. The density wave is assumed to be ordered 
{a < 0) but strongly pinned, so that it can be reduced to 
an Ising variable, [38] a [= + for c^x, - for (f)y]. Then i?^ 
becomes 



<ij> i 



(6) 



where each at represents the average on a patch, the 
first sum is over nearest neighbors (< zj >) on a square 
lattice of patches, and hi is a random variable in the 
range (— /io,/io). The strain remains a continuous vari- 
able, but now defined on patches. Since only the devia- 
toric strain rj couples to the stripes, the bulk dilational 
(ei) and shear (62) strains can be eliminated from the 
problem. However, there is a compatibility condition 
relating the strain components, to satisfy St. Venant's 
principle, which leads to a long-range interaction of the 
?7 strains, [37, 39] 



s-^ / cos [4(%)] 



where we have used the fact that ct^ is an Ising variable, 
used a renormalized 5, and 9ij is the angle between grains 
i and j measured from the Cu-0 bond direction. Numer- 
ical results for the DW ai and strain rji fields are plotted 
in Fig. 2, where we have taken (in units where a = 1) 
ft-o = 0.7, //-D^ = 8 [D is the average patch size], J = 
[since the results are not sensitive to small values of this 
parameter] , and let i5 vary. The small value of a suggests 
that the strain field is nearly unstable. This would be 
expected, since this strain couples strongly to splitting 
of the VHS peak. 

Figure 2 shows how the random patches of DW order 
generate strain fields, and how the strains develop long- 
range coherence as the DW-strain coupling is increased. 




FIG. 2: (Color online.) Co-evolution of DW and shear 
[or nematic] orders in Bi2201. Left side = maps of 
100x100 DW domains, with black dots corresponding to x- 
directed stripes {oi = +1) and blue dots to y-directed stripes 
((Ti = —1). Right side = corresponding strain fields 77, with 
magnitude given by color bars. The different rows correspond 
to 5 = 0.2 (a,b), 0.94 (c,d), 0.945 (e,f), and 1.0 (g,h). Calcula- 
tions assume periodic boundary conditions, which may limit 
domain size in frames e-h.Add system dimensions to caption 
and a, rj to the colorbox. 



Frames (a), (c), (e), and (g) show the DW domains for- 
four increasing values of S, while frames (b), (d), (f), and 
(h) show the corresponding strain fields. As i5 — >■ 0, the 
strain vanishes, but even at the low coupling S = 0.2 
in Fig. 2(b), the strain correlation length is much larger 
than the DW coupling. As 5 is increased to 0.94, the pat- 
tern of DW domains hardly changes, while the strains are 
now correlated over essentially the entire field of view, 
Fig. 2(d). For a small additional increase in coupling, 
5 = 0.945, the DW domains develop significant correla- 
tion. Fig. 2(e), and for shghtly stronger coupling both 
DWs and strain are fully ordered. In the stronger cou- 
pling regime, Figs. 2(e-h), there is a clear excess of one 
sign of DW in a given strain domain, an effect which 
has not been observed in STM[40]. In contrast, such an 
excess is weak in the intermediate coupling regime, cor- 
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responding to Figs. 2(a-d), which shows a clear similarity 
to the STM spectra of Ref. 3. 

While stripes have clearly been seen in LSCO, their 
presence in other cupratcs, inchiding the Bi cupratcs and 
CCOC, is less clear. Within our model the CDW should 
be a stable phase. In earlier STM studies the smectic 
modulations, or 'electronic glass phase', appeared clearly 
only ~30-40 meV away from the Fermi level. The modu- 
lations were correlated with the pseudogap energy scale, 
and appeared strongest at energies above 0.6 of the pseu- 
dogap energy [9], although they remain visible down to at 
least 0.2 of the pseudogap, with an energy-independent 
correlation length [3]. At lower energies a similar modula- 
tion is present, but lacking contrast reversal. [11, 41]. For 
this feature, the first harmonic ~l/4 component is actu- 
ally found to be strongest near the Fermi level. In con- 
trast, the third harmonic ~3/4 component is strongest 
near the pseudogap energy. [42] It may be that supercon- 
ductivity plays a role in making these low-energy features 
electron-hole symmetric, an effect missing from our nor- 
mal state calculation. 

While the smectic phase in Bi2212 has not been ob- 
served in x-ray diffraction or coherent x-ray scattering 
studies [43], this could be an effect of disorder, since the 
correlation length is about an order of magnitude smaller 
in Bi2212 than in LSCO. [10] Alternatively, it could be 
that strain couples the Pomeranchuk mode to a lower- 
lying mode, and level repulsion prevents the Pomer- 
anchuk wave from going soft. Interestingly, octahedral 
tilt and oxygen dimpling modes also couple to strain, 
and these modes do have low-temperature instabilities in 
many cuprates[32]. 

If the strains are confined to the Cu02 planes, then 
they would lead to a local misregistry between Cu and 
Bi atoms. A recent STM experiment has used Zn atoms 
to probe this, and the results they find are of the form 
expected for our model[44]. However, these distortions 
may be instrumental as they are not reproduced when the 
direction of the scan is changed. Alternatively, the strains 
may couple between planes, with little misregistry. We 
hope our work will stimulate further experimental work 
to determine the real strains. 

We have shown that ANN stripes provide a good model 
for the pseudogap phase in most hole doped cuprates. 
Specifically, the ANN phase in the Bi cuprates has the 
observed doping dependence of incommensurability, the 
CO displays the contrast reversal seen in STM, and the 
FS has the non-(7r/2,7r/2)-centered nodal pockets seen 
in ARPES. A secondary magnetic order enhances its re- 
semblance to a conventional stripe phase. The phase has 
quadratic coupling to a g = shear mode, which may 
be the anomalous ncmatic phase seen in STM and neu- 
tron scattering. The important role of strain in high-T 
superconducting cuprates has been noted previously. [45]. 
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A. Calculational Details 



Our investigations are based on the following hamiltonian 



H — Hf, + Hei-ph + 



(SI) 



where H^, denotes the Hubbard model, Hf,-ph the coupling between electrons and phonons 
and Hpfi the bare phonon part. In the Hubbard model 



q j destroys (creates) an electron on lattice site Ri and rij^o- = q^^Q^o-- We incorporate band 
structure effects by using for the hopping parameters tij a one-band tight-binding fit to the 
LDA dispersion in Bi-2201.^^ Interaction effects are incorporated via the (time-dependent) 
Gutzwiller approximation [(TD)GA], leading to the GA-I-RPA charge susceptibility.^^ 

With regard to the electron-phonon coupling Hei-ph we use a generic phonon model 
consisting of only longitudinal and [in-plane] transverse acoustic branches, atomic mass M, 
and electron-phonon coupling 8t. The key ingredient is that the phonons modulate the 
hopping parameter, with 'longitudinal' meaning the modulation 5t varies along the phonon 
propagation direction, and 'transverse' meaning the modulation is at right angles to the 
propagation. The corresponding operator reads as 



where 5tij = Ujaij and denotes the displacement of the atom at site Ri in direction ji. 
Finally, the phonon part is given by 



which can be diagonalized to yield the bare phonon frequencies [f^^^]^ = 2{K/M)^{2 — 
cos^qxO) — cos{qya)). Here /T^, denote effective spring constant and atomic mass for the 
longitudinal and in-plane transverse — h [T]) acoustic mode, respectively. 

In reality, a strong modulation of 5t can be produced by several phonons, including 
some involving motion of oxygen atoms perpendicular to the Cu02 planes^^. Hence in our 
model, the bare acoustic frequencies are adjusted to approximate oxygen modes in undoped 






a/3q 



aq 
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La2Cu04 [as a generic single-layer cuprate],^^ which gives {K/M)la — '^{K/M)ta = ^l^, 
taking Q° — 12.4 meV, and M is the oxygen mass. In this case, the dominant g'-vectors, 
Fig. SI, should be approximately correct, since they are controlled by Fermi surface nesting, 
but the 'phase diagram'. Fig. SI, is merely indicative of regions of strong clcctron-phonon 
coupling and significant Kohn anomalies. A realistic calculation of instability would require 
a model of the full phonon dispersion, and in particular, how the strong softening of one 
branch interacts with other branches of the same symmetry at lower energy. 

The electron-phonon coupling leads to dressed phonon frequencies according to 



(S2) 



where the renormalized elastic constants can be evaluated within the TDGA as^^ 



(S3) 



In terms of a bare susceptibility Xoq = Xok,q, with 



The terms in Eq. S3 are 
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with £'k,q,(T = ^k+q,CT + ^k,<T' ^^^^ dispcrsiou e^, dressed dispersion ek 
Gutzwiller renormalization factor Z — z^. Finally, 



(S7) 

Ze^, with the 



Jk,k+q,iJ, ~~ ^'''^0Jk,k+q,ix^ 



2 AO) 



(S8) 



fkl+q,iJ.Qq is Fourier transform of the SSH interaction /j^,^ = tijaij{u'^ - uf), fq^iQ'^ 
is the Fourier transform of fi^^ — fij,^, Qg is the Fourier transform of uf, aij = 
— din (tij)/dr > and fi = x,y. Equation S2 can be written in terms of a Stoner fac- 
tor, which takes a simple form when 5K is diagonal [as along F ^ (tt, 0), F ^ (tt, tt)], 

UeffXo^ = -5K,jM,[nl]\ (S9) 

It is convenient to normalize = MnE^SK^^/j,, so that the strength of the electron- 

phonon coupling is measured in terms of 

where for nearest neighbor hopping, ctjj — ^/a and a is the Cu-Cu separation, and we have 
introduced the electron-phonon coupling constant 

Taking M to be an oxygen mass, E*j^o = Er/^^ = 16.7 pieV. For Bi2201, with t = 435 meV, 
this is equivalent to Agpo = — 0.047. Note that Agp is independent of the phonon 

mass. 



B. Results: Soft phonons, Charge stripes, and Phase Diagrams 

Figure SI compares the bare and renormalized LA and [in-plane] TA phonon frequencies 
in Bi2201. Along the (tt, 0) (tt, 7r)-branch the modes are mixed, and labelled as predom- 
inantly longitudinal or transverse. The sharp dips in the dressed frequencies are caused 
by peaks in the bare susceptibihty associated with FS nesting. Each peak in x leads to a 
prominent Kohn anomaly in the phonon spectrum, which can lead to an instability if the 
renormalized fi^^j becomes negative. By comparing the present results with earlier calcula- 
tions for magnetic stripes, we find that the instabilities fall at nearly the same g-values for 
both kinds of stripes as a function of doping, being controlled by the same Fermi surface 
nesting. 

In the doping range appropriate to the cuprates, there are two competing instabilities. 
For very low doping the dominant instability is to an incommensurate state at the edge of 
the (tt, tt) plateau at (tt — 5, tt) [anomaly 2 in Fig. SI (a) - called a plateau instability], while at 
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higher dopings a diagonal ANN phase is dominant at (tt — 5, tt — 5) [4 in Fig. Sl(a)]. Figure SI 
compares the phonon softening in a weak couphng {U/Ubr — 0.2) and an intermediate 
couphng case (U/Ubr = 0.6) more appropriate for cuprates. 

Figure S2(a) shows the resulting charge order phase diagram for several values of U. As 
expected, when U is large (> Ubr, the Brinkman-Rice coupling) charge order is strongly 
suppressed at half filling, but it is restored rather quickly when U is reduced or x increased. 
ANN stripes (circles and squares) are dominant for a wider doping range than in the mag- 
netic phase diagram, but near half filling and in the overdoped regime vertical/horizontal 
stripes associated with the (tt, 7r)-plateau (triangles and diamonds) win out. Near the VHS 
{x = 0.42) the value of the critical parameter goes to zero, and it is found that the domi- 
nant nesting g-vector approaches F, corresponding to phase separation. There is an overall 
similarity with the magnetic-order phase diagram. Fig. S2(b), found earlier. Similar phase 
diagrams have been found for most cuprates studied, except for the magnetic phase diagram 
of LSCO, for which the VHS is much closer to half filling, resulting in the dominance of the 
(tt, 7r)-plateau instabilities. 

Figures 1 and S2 give a very graphic picture of the role of Fermi surface nesting in the 
cuprates: not only are the dominant g- vectors correctly predicted, but also the correct phase 
between two competing DW instabilities. However, exact agreement between the FS nesting 
vector and the experimental DW vector is not to be expected, even for conventional CDWs 
and SDWs. First, nesting is a manifestation of a Stoner-type instability, l — Uxoi'i,^ = 0) = 
0, where U is here a Gutzwiller-modified Hubbard U and Xo is a bare [charge or magnetic] 
susceptibility. Only at threshold U = Uc is FS nesting exact. For U > Uc the dominant 
g- vector adjusts to further lower the electronic energy. Numerical estimates suggest the shift 
is relatively small unless the FS is fully gapped. One common effect in conventional DWs is 
a commensur ability pinning: the DW can gain energy by shifting from an incommensurate 
nesting vector to a nearby commensurate vector. This is likely to be important in the 
cuprates, where the DW periodicities correspond to only a few atoms across. For instance, 
several cuprates seem to have the DWs pinned near a 4-Cu periodicity, and the Bi-2201 data 
in Fig. 1 resemble a crossover from period 4 to period 6 [near nesting] with increasing doping. 
For Bi-2212, the g'-vector is more consistent with the prediction for the bonding band. It 
may be that proximity to period-4 order stabilizes the bonding band DW over the competing 
antibonding band DW. Another factor is that in two-dimensions the lowest-energy DW state 
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corresponds to double nesting - simultaneously nesting two disconnected FS sections. This 
can be achieved either by a single DW or by two or more coherent DWs^^. As noted above, 
the latter seems to be the case for the ANN phase in the Bi-cuprates: the spatial pattern 
seen in Bi-2201^ can be fit as a coherent superposition of horizontal and vertical stripes. 
The superposition of orthogonal CDWs forms a distinct pattern from a checkerboard. See 
Fig. 1(c) of Ref. [6]. Finally, even when the DW represents a local minimum of the free 
energy, the global minimum may correspond to a state of [nanoscale] phase separation, as 
predicted in Bi-2212.^2 

C. Period 6 Stripes 

Further insight into the ANN phase comes from comparisons with photoemission and 
tunneling data. This requires real space modeling, as in Fig. S3. We adjust the doping to give 
a simple commensurate super lattice vector Qq. For definiteness we model Bi2212, neglecting 
bilaycr splitting (dotted line in Fig. 1). In this case Qq is exactly 27r/6a, where a is the lattice 
constant, at a doping x = 0.19. We model the onc-dimcnsional ANN order as a Pomeranchuk 
wave, with modulated hopping parameter tx — t ± St, ty — t with St proportional to the 
lattice distortion and the maximum St — O.lt, neglecting possible modulations of t', t", and 
t'". We assume a bond-centered lattice distortion. Fig. S3(c). Since the STM images look 
more hke ID stripes in Bi2212, and crossed (2D) stripes in Bi2201, we model both, with 
the ID results in Fig. S3(a,d,e) and the 2D results in Fig. S3(b,f,g). For this distortion. 
Figs. S3(a,b) show the density-of-states (DOS), both the average value (red lines) and the 
local values on different sites. There is a local pseudogap on the stripes, of a magnitude 
comparable to that seen in the STM experiment. For the ID stripes. Fig. S3(a) shows that 
the short-bond rows have a larger local DOS for filled states. There is a clear contrast reversal 
between electron and hole states - a higher DOS for filled [empty] states on the large-t [small- 
t] atoms. Similar results arise for the 2D patterns, Fig. S3(b). An approximate contrast 
reversal is commonly found in CDW systems^^, and the resulting pattern bears a striking 
resemblance to the phase seen in STM experiments. ^'^'^^'^^ Figures S3(d,f ) show the resulting 
FSs, with a structure factor (SF) correction appropriate for angle resolved photoemission 
(ARPES) spectra^^'^. The spectral weight is plotted on logarithmic scale to enhance weak 
superlattice features. When the SF is included, it is seen that most of the spectral weight 
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lies near the original FS. For either ID or 2D stripes, the FSs are composed of nodal arcs, 
as seen in experiment, with weaker features near the antinodes. 

The present calculation represents a minimal model of the stripes. We note that a number 
of refinements are possible. First, in principle one could also include an on-site (Madelung) 
correction of ~+[-]50 meV per short [long] bond, as holes are attracted by the negatively 
charged oxygens. In the absence of this term, we find that holes tend to pile up on the 
atoms with more long bonds, and avoid the atoms with short bonds. Moreover, a number 
of secondary order parameters are not only possible, but likely, greatly complicating the 
DW picture. As can be seen from the FS maps, Figs. S3(d,f), there are a number of level 
crossings that the Pomeranchuk wave has not gapped. Additional electronic distortions can 
take advantage of this, opening gaps at the crossing points to lower their energy and further 
stabilize the distortion. 

Here we give a single illustration of this effect. While this is a conventional charge density 
wave (CDW), in the cuprates it can acquire aspects associated with stripe physics. Thus, 
rows with lower doping will display enhanced Mott physics while the hole-doped rows would 
be closer to optimal doping. In Fig. S3(c) we show an example where proximity to half filling 
induces a secondary AFM order on the less doped stripes. Fig. S3(e) shows the resulting FS 
map in the ID case. Remarkably, the Fermi arcs have become pockets (arrow), which are 
not centered on (7r/2, 7r/2), as observed in some experiments. Similar effects arise in 2D, 
Fig. S3(g), but here optimal doping shifts to higher x — 0.31 to better nest the antinodal 
regions. Note that the Fermi arc seems to terminate along the AF zone boundary, consistent 
with STM^, while the FS contour is quite squarish, which has been taken as a signature of 
underlying stripe order in ARPES studies. 

Recent quantum oscillation (QO) studies have found evidence for small FS pockets in both 
hole-doped YBa2Cu307-5 (YBCO)Sio-i2 and electron-doped Nd2-^Ce^Cu04 (NCCO)^^^ 
While the present model successfully predicted the areas in NCCO^^^, earlier proposed 
pggSi5,i6 YBCO bore little resemblance to the present results. However, after this 
manuscript was completed, we became aware of new results which found that there is most 
likely only a single (electron-like) pocket in underdoped YBCO^^^, which can be modeled 
by a two-q density wave^^^, similar to the present result. While the doping range in YBCO 
where these QOs are observed probably corresponds to period-4 stripes, the proposed Fermi 
surface preserves the nodal states, and overall bears a close resemblance to the high spectral 
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weight star-shaped feature in Fig. S3f. Indeed, we find that a similar nodal Fermi surface 
exists in Bi-2212 in this period-6 crossed stripe state. 

Finally, we note that ANN nesting provides a natural explanation for one more feature 
of the ARPES data. Yang, et al.^'' noted that the pseudogap in Bi2212 is centered at the 
Fermi level in the antinodal regime. As they point out, this is unexpected for a simple 
(tt, tt) antiferromagnetic order. However, since the ANN phase specifically involves nesting 
along the long fiat sections of the FS in the antinodal regime, it is natural to expect that 
this nesting is centered at the Fermi level. We show that this is in fact the case in Fig. S4. 
Figure S4(a) shows the dispersion along a cut in the antinodal regime, corresponding to the 
light blue line in Fig. S3(f). It can be seen that the ANN gap is centered on the Fermi 
level. Further confirmation is shown in Fig. S4(b), which plots the spectral weight along the 
dotted fine in Fig. S4(a). To compare with experiment [Fig. 4(f) of Ref. [17]], the data in 
this figure have been broadened by 7 = 0.33a; + 5 meV. Clearly the ANN nesting provides 
a natural (non-superconducting) explanation for this feature. 

Half pockets, as in Fig. S3(e), are also found in a phenomenological model of the cuprates 
based on the t — J model. However, this model cannot explain CO with nesting vectors 
different from (tt, tt), and must explain the pseudogap centered at the Fermi level as a 
consequence of a superconducting gap.^^ But this is inconsistent with growing evidence 
that superconductivity only appears at an onset temperature T^^get < T*, the pseudogap 
temperature. ^^'^"^^ 

The present model bears some resemblance to the valence bond glass (VBG)^^^'^^. Indeed, 
since impurities are a relevant perturbation of a CO phase, the ANN phase will become 
more glass-like in a real Bi-2212 sample with large interstitial-oxygen disorder. However, a 
significant difference is that the one-dimensional ANN phase has no bond length modulation 
along y (the direction along the stripe), while the VBG has a d- wave-like modulation. 

I. SUMMARY OF MODEL CALCULATIONS 

The present letter plus supplementary material describes three closely related calculations 
which together provide a coherent microscopic model for charged stripes, to complement our 
earlier work on magnetic stripes. The three calculations are designed to answer three ques- 
tions: (1) what is the phase diagram for predominantly charge-based stripes? (2) can these 
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stripes explain the various features of stripes seen in ARPES, STM, and quantum oscillation 
measurements? and (3) how are the 'nematic' phenomena related to these stripes? Each 
calculation requires progressively more approximations, so in this subsection we summarise 
them. 

The underlying microscopic model contains three key ingredients. (1) The dominant 
electron-phonon interaction arises from modulations of Cu-Cu bondlength, which modify the 
electronic hopping energies. (2) Since the energies change as high powers of the bondlengths, 
there is also a second-order (t^) coupling to the electronic nematic order, which couples 
directly to elastic shear strain. (3) The strain coupling is long range, so that even if charge 
pinning restricts the stripes to nanoscopic length scales, the strain-nematic order can persist 
over substantially larger length scales, as found in experiment. We have amply documented 
all three ingredients, so that the model qualitatively describes the observed phenomena. The 
question remains, to what extent is the agreement quantitative? 

For the phase diagram calculation [Figs. 1, SI, and S2], the electronic parameters [hopping 
energies and Hubbard U] are the same as used for our magnetic calculations, and should 
be quite accurate. The electron-phonon interaction parameter 7 is consistent with litera- 
ture values, although it is not clear how the 7s vary for more distant hoppings, so we have 
simply neglected higher order terms. The main assumption is in the bare phonon frequency 
Jlo(?)- We have introduced an effective one-band model, to approximately describe the Cu 
displacements and corresponding in-plane 'acoustic' modes. A more realistic model would 
be desirable, but would require a model of the full 'bare' phonon dispersion, including eigen- 
functions for all branches, in the absence of this electron-phonon coupling term. We believe 
that the present model properly describes the competition between Hubbard-t/ suppression 
of charge order and the tendency of electron-phonon coupling to form stripes, and captures 
the dominant Fermi-surface nesting vectors. The approximations on the phonon spectrum 
will mainly affect the exact coupling 7 for instability and the nature of the phonon mode that 
condenses. For example, in LSCO the low-temperature orthorhombic and low-temperature 
tetragonal phases both involve soft modes of octahedral tilts. Interestingly, both modes 
couple to shear strains, and both do go unstable in appropriate doping ranges. 

The phase diagram gives the minimum 7 (or equivalently, electron-phonon coupling A) 
needed to produce a CDW order, and the corresponding g'-vector. To find the corresponding 
real and k-space stripe structure [Figs. S3, 4] we should choose a 7 larger than the threshold 
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value and numerically minimize the free energy on a lattice, which will produce a shifted q- 
vector, as found in the magnetic stripe case, Ref. [16]. Instead, we simply assumed a plausible 
amplitude for the lattice distortion and calculated the corresponding electron energy levels. 
We also estimated the effect of a secondary magnetic order, as expected for stripes near half 
filling. We believe that such approximations are adequate to capture the qualitative features 
of the experiments. 

Finally, the nematic order arises naturally in the presence of a shear strain which splits 
the VHS degeneracy. Wc demonstrated how a shear strain arises as a secondary order 
parameter in the presence of stripes. The emergent nature of the nematic phase arises 
from the long-range nature of the strain coupling, which makes it much more robust than 
the stripe order in the presence of point impurity disorder. In calculating Fig. 2 we were 
interested in one issue: could the nematic order have a much larger correlation length that 
the underlying stripe order, in the presence of impurity pinning? In the main text we have 
provided a concrete calculation showing that the answer is yes, while above we presented a 
possible microscopic realization. 
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FIG. SI: Bare phonon dispersion (dashed line) compared to dressed dispersion assuming 
U/Ubr = 0.20 (light lines) or 0.60 (dark lines) at a series of hole dopings x = (a) 0.05, (b) 
0.10, (c) 0.20, (d) 0.30. Longitudinal [transverse] phonons in shades of blue [red]. Mate- 
rial parameters appropriate for Bi-2201, for which Ubr = 13. 6t. Only modulation of the 
nearest-neighbor hopping t is included, with doping independent magnitude 7 = 3.28. By 
convention, real fip^'s are plotted as positive numbers, imaginary fip^'s as negative. 
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FIG. S2: Phase diagrams for Bi2201. (a) Charge order, showing threshold strength of 
electron-phonon coupling Agp for several values of Hubbard U. Symbols represent different 
stripe symmetries, as vertical [diagonal] (vr, 7r)-plateau stripes = triangles [diamonds]; ver- 
tical [diagonal] ANN stripes = squares [circles]. Dashed lines indicate transitions between 
different symmetries, while dotted line corresponds to 7 = 3.5. (b) Magnetic order phase 
diagram, from Ref. 1. 
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FIG. S3: ANN superlattice for Bi-2212. (a,b) Average (red) and local dos for ID (a) and 
2D (b) stripes. Defining s and / for short and long bonds, a site with two short bonds as 
(s, s), etc., the various local dos are represented as follows: in (a): {s,s) (blue), (s, /) (light 
green), and (/, /) (green). In (b): (s, s)^, (s, s)y (dark blue), (s, s)^, (s, l)y (blue), (s, s)^, (/, l)y 
(light blue), {s,l)^,{s,l)y (orange), {s,l)^, {l,l)y (light green), and (1,1)^, {l,l)y (green), (c) 
Local distortion pattern for ID stripes. Also shown is pattern of AFM order assumed for 
Figs. S3(f,g). (d,f) Fermi surface maps for ID (d) and 2D (f) stripes, with spectral weight 
plotted on a logarithmic scale. (e,g) Corresp^(^ding FS maps with added AFM order. Arrow 
in (e) illustrates a prominent half-pocket Fermi surface section. 
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FIG. S4: {a) E — k map of electronic dispersion along the cut in Fig. S3(f) [blue line], 
showing ANN gap centered at the Fermi level Ep. (b) Plot of spectral weight vs binding 
energy along the dotted line in frame (a). For (b), an energy dependent broadening 
7 = 0.33a; + 5 meV was included to mimic the experimental results. 
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